Pairing fluctuation effects in a strongly coupled color superfluid/superconductor 
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We investigate the effects of pairing fluctuations in fermionic superfluids/superconductors where pairing oc- 
curs among three species (colors) of fermions. Such color superfluids/superconductors can be realized in three- 
component atomic Fermi gases and in dense quark matter. The superfluidity/superconductivity is characterized 
by a three-component order parameter which denotes the pairing among the three colors of fermions. Due to the 
SU(3) symmetry of the Hamiltonian, one color does not participate pairing. This branch of fermionic excitation 
is gapless in the naive BCS mean-field description. In this paper, we establish a pairing-fluctuation theory for a 
strongly coupled color superfluid/superconductor and investigate the pairing fluctuation effects on the unpaired 
color. At low temperature, a large pairing gap of the paired colors suppresses the pairing fluctuation effects 
for the unpaired color, and the spectral density of the unpaired color shows a single Fermi-liquid peak, which 
indicates the naive mean-field picture remains valid. As the temperature is increased, the spectral density of the 
unpaired color generally exhibits a three-peak structure: the Fermi-liquid peak remains but get suppressed, and 
two pseudogap-like peaks appears. At and above the superfluid transition temperature, the Fermi-liquid peak 
disappears completely and the all three colors exhibits pseudogap-like spectral density. The coexistence of Fermi 
liquid and pseudogap behavior is generic for both atomic color superfluids and quark color superconductors. 
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I. INTRODUCTION 



Color superconductivity in dense quark matter in general 
appears due to the attractive interactions in certain diquark 
channels OJ-dl] . Taking into account only the screened (color) 
electric interaction which is weakened at the Debye mass scale 
g/j. (g is the QCD coupling constant and ji is the quark chem- 
ical potential), the early studies |Q]] predicted a rather small 
pairing gap A ~ lMeV for moderate density where fi ~ Aqcd 
(Aqcd ~ 300MeV is the QCD energy scale). The break- 
through in this field of research was made in J2l where it 
was observed that the pairing gap is about 2 orders of mag- 
nitude larger than the previous prediction using the instanton- 
induced interactions and the phenomenological four-fermion 
interactions. On the other hand, it was first pointed out by 
Son that at asymptotic high densities the unscreened mag- 
netic interaction dominates 0]. This leads to a non-BCS gap 
A ~ /ig~ 5 exp (-c/g) JH where c = 3n 2 / V2, which matches 
the large magnitude of A at moderate density predicted by 
the instanton-induced interactions and the phenomenological 
four-fermion interactions. 

Such gaps at moderate density are so large that they may 
fall outside of applicability range of the usual BCS-like mean 
field theory. It was estimated that the size of the Cooper pairs 
or the superconducting coherence length f c becomes com- 
parable to the averaged inter-quark distance d ^ at moder- 
ate density where yt ~ Aqcd- This feature is highly con- 
trasted to the standard BCS superconductivity in metals where 
f c » d. Qualitatively, we can examine the ratio of the su- 
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perconducting transition temperature T c to the Fermi energy 
E F , k = T C /E F 01- We have k ~ 1(T 5 for ordinary BCS 
superconductors and k ~ 1CT 2 for high temperature super- 
conductors. For quark matter at moderate density, taking 
E F a 400MeV and T c 50MeV H, we find that k is even 
higher, k ~ 10 _1 , which is close to that for the resonant super- 
fluidity in strongly interacting atomic Fermi gases 071]. This 
indicates that the color superconducting quark matter at mod- 
erate density is in the strongly coupled region or the BCS- 
BEC crossover JHOiil]. It is known that the pairing fluctuation 
effects play important role in the BCS-BEC crossover lfl2tl . 

Many efforts have been made in understanding the single- 
particle properties and equations of state in the strongly inter- 
acting region of the s-wave Fermi superfluids lfl3l \\M where 
pairing occurs between two components of fermions. Satis- 
factory agreement with the experimental data from trapped 
fermionic atoms has been achieved lfl5l [l6ll . Recently, the 
single-particle spectral density in the strongly interacting re- 
gion has been measured in cold atom experiments. It was 
found that a pseudogap phase appears above the critical tem- 
perature lfl7ill8ll . where the system retains some of the charac- 
teristics of the superfluid phase such as a BCS-like dispersion 
and a partially gapped density of states but does not exhibit 
superfluidity. The pseudogap behavior of the single-particle 
excitations was also found by Quantum Monte Carlo calcu- 
lations 11911 . The agreement between the theoretical predic- 
tions and the experimental data in the equations of state and 
the single-particle spectral functions indicates that the pairing 
fluctuation effect is significant in the strongly interacting re- 
gion, which defies the conventional BCS-like mean field the- 
ory. 

In this paper, we investigate the pairing fluctuation ef- 
fects in color superfluids/superconductors where pairing oc- 
curs among three species (colors) of fermions, in contrast to 
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three-component Fermi gas and (ii) a two-flavor quark color 
superconductor at quark chemical potential jj ~ 400MeV. 

The beyond-mean-field approach for the study of the pair- 
ing fluctuation effects adopted in this paper is a generalization 
of the T-matrix theory for two-component (spin-j) fermionic 
systems developed by Pieri, Pisani and Strinati [13]. We 
will focus on the pairing fluctuation effects on the unpaired 
fermions in the color superfluids/superconductors, which is 
absent in the two-component systems. The paper is organized 
as follows. In Sec. II, we study the pairing fluctuation ef- 
fects in atomic color superfluids, which may be realized in 
three-component atomic Fermi gases. In Sec. Ill, the pair- 
ing fluctuation effects in two-flavor color superconductors will 
be discussed. We summarize in Sec. IV. The natural units 
c — H — ks = 1 will be used throughout the paper. 



FIG. 1 : (color-online) A schematic plot of the pairing pattern in color 
superfluids/superconductors, corresponding to the order parameter 
gauge A = (A, 0, 0). Only the red and green fermions participate 
the pairing and their pairs condense. Even though the blue fermions 
do not participate pairing, pairing fluctuations exist for red-blue and 
green-blue pairs, in addition to the red-green pairs. 



the ordinary BCS-BEC crossover problem where pairing oc- 
curs between two fermion components. Such color super- 
fluidity/superconductivity can be realized in three-component 
atomic Fermi gases Il20ll and two-flavor dense quark matter 
which appears around ju ~ 400MeV in the quark matter phase 
diagram 12111 . 

The color superfluidity/superconductivity is generally char- 
acterized a three-component order parameter A = (Ai,A2, A3), 
where Ai ~ (iff T if/ g ), A2 ~ {iffgftb), an d A3 ~ (i/z^b)- Here 
we use red (r), green (g), and blue (b) to denote the three 
fermion species. In QCD, they are the color degrees of free- 
dom for quarks. If the Hamiltonian of the system has a SU(3) 
color symmetry, the non-vanishing order parameter breaks the 
SU(3) symmetry group down to a subgroup SU(2). However, 
due to the SU(3) symmetry of the Hamiltonian, the effec- 
tive potential *V(A) should depend only on the combination 
AA' = |Ai[ 2 + IA2P + IA3I 2 . Therefore, we can choose a spe- 
cific gauge for the three-component order parameter, such as 
A = (A, 0, 0) without loss of generality. This simple argument 
shows that there exists a branch of fermion which does not 
participate the pairing. For the gauge A = (A, 0, 0), only the 
red and green fermions participate the pairing, while the blue 
fermions do not. A schematic plot of this pairing pattern is 
shown in Fig. Q] In the naive BCS-like mean-field descrip- 
tion, the unpaired blue fermions possess a free energy disper- 
sion and is gapless. However, we note that even though the 
red-blue and green-blue pairs do not condense, their pairing 
fluctuations do exist. For weak attraction, the pairing fluctua- 
tion effects on the single-particle excitations can be safely ne- 
glected and the BCS-like mean-field description is applicable. 
In this paper, we are interested in the strongly interacting sys- 
tems, where the pairing fluctuation effects become significant. 
Two systems will be considered: (i) a resonantly interacting 



II. COLOR SUPERFLUIDITY IN ATOMIC FERMI GASES 

In this section we study the pairing fluctuation effects in 
an atomic color superfluid B22H28I1 . a cold atom analogue of 
color superconductivity in dense QCD l29ll . We consider a 
dilute Fermi gas composed of three species of fermions with 
a common mass m. Such a system can be realized in cold 
atom experiments by trapping the lowest thre e hy perfine states 
of the 6 Li or 40 K atoms in the atomic trap 12011 . We assume 
that there exist short-range attractive interactions among dif- 
ference species. The attraction strength can be tuned from 
weak to strong by means of the Feshbach resonance. In the 
dilute limit, the attractive interactions can be modeled by con- 
tact interactions. The Hamiltonian density of the system is 
given by 



S>a=l 



where ^1,^2, and //3 are the chemical potentials for the three 
species, and g\2,gn,gn are the contact interactions among 
them. 

In the following, we assume that the total particle number 
is fixed and the chemical potentials become equal, fii — fj.2 — 
yt/3 = ju. Further, we assume that the three coupling constants 
also equal each other, g\i = gn = g23 = g- Then the Hamilto- 
nian (Q]) has a global SU(3) symmetry. To show this explicitly, 
we define a three-component fermion field iff = , 1J/7, i^3) T . 
Then the Hamiltonian density ([1) can be expressed as 112411 

^ = ^ " ^ * + 4 2 (^a^X^alA), (2) 

where A d are the Gell-Mann matrices in the "color space" 
spanned by the three fermion species. The contact coupling 
constant g can be renormalized by introducing the s-wave 
scattering length a s of the short-range potential via the equa- 
tion 
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Here = k 2 /(2m) is the free dispersion of the fermions, and 
the integral over the momentum k is associated with a cutoff 
A. We can set A — > oo in the physical equations since all UV 
divergences are removed by Eq. 0). In general, by tuning the 
scattering length from small to large negative values, we can 
realize an evolution from a weakly coupled BCS-like color 
superfluid to a strongly coupled color superfluid. 



A. BCS mean-field theory 

Due to the attractive interactions among the unlike col- 
ors, at low temperature the system is in a color superfluid 
phase. Different to the two-component (spin-^-) fermionic 
systems where the superfluidity is characterized by a single- 
component order parameter A ~ (ifrfif/i), the color superflu- 
idity in the present three-component Fermi gas is character- 
ized by a three-component order parameter A = (Ai, A2, A3), 
where 

A 2 = !<<A T /i 5 <A> = fe<toto>. 



A 3 = |<-A T i 7 <A) = /g<<A3<Ai>. 



(4) 



Since the Hamiltonian ([2]) has an exact SU(3) symmetry, we 
can show that the effective potential 'V(A) depends only on 
the combination [24] 



(5) 



AA^ = I Ail 2 + IA2I 2 + I A 3 1 2 . 



Therefore, different pairing configurations are physically 
equivalent and we can choose a specific gauge A t - A ^ 
0, A2 = A3 = without loss of generality. In this gauge, 
only the red and green fermions participate the pairing and the 
red-blue pairs condense, leaving the blue fermions unpaired. 
For a general gauge where all three components are nonzero, 
the unpaired branch is a linear combination of the three col- 
ors. Due to the SU(3) symmetry of the hamiltonian, different 
choices of the order parameter lead to the same physical re- 
sults. For convenience, we use the gauge A = (A, 0, 0) in the 
following. 

Before going on, we should mention that the BCS-BEC 
crossover in the present three-component Fermi gases is 
somewhat different from the two-component systems due to 
the possibility of the appearance of a trionic phase, which is a 
Fermi liquid state of three-fermion bound states 12311 . Based 
on an attractive Fermi Hubbard model, it was shown that there 
exists a quantum phase transition from the color superfluid 
phase to the trionic phase when the attraction strength exceeds 
a critical value l23ll . On the other hand, it was shown that large 
three-body losses in three-component Fermi gases confined in 
optical lattices can stabilize the color superfluid phase by sup- 
pressing the formation of trions l27ll . Therefore, in this paper 
we neglect the possibility of the trionic phase and consider the 
color superfluid phase only. 

Now we turn to the general formalism. It is convenient to 
work with the Nambu-Gor'kov basis ^ng = (ft, | A*) T - In the 



Nambu-Gor'kov representation, the fermion self-energy 2,(K) 
and the dressed fermion propagator S{K) are 2 x 2 matrices. 
They satisfy the Dyson's equation 



S n (K) Sn(K) 
S 21 (K) S 22 (K) 

W>-& WZ„(Z) E 12 (£) 

iCOn+fr J I £21 (X) Z 22 (K) 



(6) 



Here and in the following, K = (ico n , k) with a>„ = {In + l)nT 
(n integer) being the fermionic Matsubara frequency, and ^ = 
ek - fi. From the Green's function relation we have the gap 
equation 



1 K 



and the number equation 

n= VTr[5ii(Z)]. 



K 



(7) 



(8) 



Here and in the following the notation Y,k - ^ Z» Zk with 
Ek = / d 3 k/(27r) 3 is used. 

In the BCS mean-field theory, the fermion self-energy is 
chosen as 



UK) = £ BCS (X) = 



AA 2 
AA 7 



(9) 



and hence momentum-independent. Here we have set A to 
be real without loss of generality. Then the dressed fermion 
propagator S{K) can be evaluated as 

S^ S (K) = g A (K)A Ig +g Q (K)A h , 
S^ S (K) = Ta(K)A 2 , 



S^ S (K) 

:>BCS 



-^ CS (-^), 
BCS/ 



S^(K) = S^(K). 



(10) 



Here we have defined two matrices A Ig = diag(l, 1,0) and 
Ab = diag(0, 0, 1) in the color space. The Green's functions 
Qh{K), Th{K) and Qq(K) are analytically given by 



0a(K) = 



(iu n ) 2 - El ' 



TdK) = - 

(ico„Y - Er 



1 



ia> n - £ k 



(11) 



where E^ = ^ 2 + A 2 is the BCS-like dispersion. 

In the above BCS mean-field description, it is clear that the 
paired colors possess BCS-like dispersions and obtain an ex- 
citation gap A, while the unpaired color has a free dispersion 
and remains gapless (for /i > 0). The physical value of the 
pairing gap A is determined by the BCS gap equation 



4na s 4^ 



1 ~ 2f(E k ) 1 
2e k 



2£ k 



(12) 
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where f(E) = l/(e E/T + 1) is the Fermi-Dirac distribution. 
However, we note that there exists a serious problem for the 
naive BCS mean-field approach. If we consider another sys- 
tem with the couplings gu = g and go = g23 = 0, i.e., the 
blue color is really free, we find no difference in comparison 
with the present system with gi2 = gn = g23 = g- There- 
fore, the pairing fluctuation effects is likely significant in the 
color superfluid, especially when the attractive coupling g is 
not weak. 



B. Fermion self-energy beyond BCS 



Q-K Q-K 



(a) 



D(Q) 



(b) L (K) 



Q-K 




(K) 






D(Q) 
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Now we take into account the pairing fluctuations and study 
their effects on the single-particle excitation spectra. To this 
end, we first construct the particle-particle ladder D(Q) or the 
"pair propagator", which is diagrammatically represented in 
Fig. |2ja). In the color superfluid phase it is a 2 x 2 matrix, 



x£(Q) xf 2 (Q) 



(13) 



where x(Q) is trie P a i f susceptibility. Note that the matrix 
elements of T>(Q) and x(Q) aie a l so matrices in an adjoint 
space of the SU(3) group spanned by the indices a,b = 2, 5, 7. 
The explicit form of the pair susceptibility x(Q) is given by 



y ab 



(0 = —-\j]TT[A,Sn(Q-K)A h S u (K)], 

Z K 

xfiiQ) = Xu(-Q), x? l (G)=Xn(Qy d4) 



Here and in the following, Q = (iv n , q) with o>„ = 2nnT (n 
integer) being the bosonic Matsubara frequency. The fermion 
self-energy beyond the BCS approximation (O can be consid- 
ered by taken into account the self-energy l. L (K) shown in Fig. 
|2jt>). Now the full fermion self-energy in our consideration is 
given by 



where ^(K) can be expressed as 



(15) 



a,b Q 



Z L l2 (K) = -^^^(0^12(2-^, 

a.b Q 

I. L 22 {K) = --L L n {-K), 5^(X) = Ef 2 (Z). (16) 

In general, the above equations together with the gap and 
number equations ^ and ([8]l form a closed set of integral 
equations for the superfluid order parameter A, the chemical 
potential fi, and the dressed fermion propagator S(K). How- 
ever, to keep the Goldstone's theorem and recover the correct 



FIG. 2: The diagrammatic representation of the particle-particle lad- 
der D(Q) and the self-energy ~L L (K). The dashed lines denote the 
coupling constant g. 



pair excitation spectrum at strong coupling, we adopt the fol- 
lowing prescriptions suggested in 11311 : (i) In evaluating the 
pair susceptibility x(Q) an d me self-energy T. L (K), we use 
the the fermion propagator S(K) of its BCS mean-field form 
([Tol l; (ii) The off-diagonal fermion propagator Sn(K) in the 
gap equation (0 is also replaced by its BCS mean-field form, 
while in the number equation ([8]l we take into account the 
pairing fluctuation effects on the diagonal fermion propaga- 
tor Sn{K). Therefore, the gap equation (0 for the superfluid 
order parameter A still takes its BCS form ( fT2l i. which ensures 
a gapless pair excitation spectrum, i.e., the Goldstone's theo- 
rem. The beyond-mean-field corrections for the pairing gap A 
and chemical potential are reflected in the full number equa- 
tion ([8]). The pairing fluctuation effects on the single-particle 
excitation spectra are included in the self-energy Z L {K). 

After some simple matrix algebras, we find that the pair 
susceptibility x(Q) an d the particle-particle ladder D(Q) are 
both diagonal in the adjoint space, i.e., 



(17) 



For the a = 2 sector which represents the red-green pairing, 
the pair susceptibility x 2 (Q) can t> e evaluated as 



*?i(0 = ~r--J\ TY&A(K)g A (Q-K)-±- 
XvSQ) = Yj T ^ k)T " {Q - K) - 



(18) 



The corresponding particle-particle ladder D 2 (Q) reads 



£>!l(0 = 



£>12(0 



xU-Q) 



X 2 n (Q)X 2 n(-Q)-\X 2 nm 2 ' 

xUQ) 



(19) 



x 2 u{Q)x 2 u(-Q)-\x\ 2 m r 

Using the BCS gap equation (fT2] i. we can show that 

^(0,0) =^(0,0). (20) 
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Therefore, the pair excitation from the a = 2 sector is 
gapless, corresponding to one broken generator A Ig . We 
note that this Goldstone mode is essentially the same as 
the Anderson-Bogoliubov mode in the conventional two- 
component fermionic superfiuids 13011 . It possesses a linear 
dispersion u>(q) oc |q| in the low momentum and frequency 
limit. In the strong coupling limit, it recovers the Bogoliubov 
excitation spectrum for a weakly interacting Bose gas 1 1 3L 1 3 Qfl - 
The sectors a = 5 and a = 7 represent the red-blue and 
green-blue pairings, respectively. They are degenerate, corre- 
sponding to the residual SU(2) symmetry group. Since these 
pairs do not condense according to our order parameter gauge, 
the off-diagonal components vanish. We have 

XniQ) = X[ 2 (Q) = 0, £>i 2 (0 = D] 2 (Q) = 0. (21) 

The non-vanishing diagonal component of the pair suscepti- 
bility can be evaluated as 



2e k 



.(22) 



The particle-particle ladder is also diagonal and is given by 

1 



2>ii(G) = 



a = 5,7. 



x d u(Q) 

Using the BCS gap equation ( fT2l , we find that 



*Ti(0,0) = 0, 



5,7, 



(23) 



(24) 



We can neglect this off-diagonal contribution since it is gener- 
ally much smaller than the BCS self-energy S BCS lfl3ll . How- 
ever, this approximation has no effect on the spectrum of the 
unpaired color we are interested in in the following. 

It is worth comparing the self-energy 2< rg (K) with that in 
two-component (spin-^) system. In two-component system, 
the self-energy T. rg (K) reads lfl3ll 



Q 



K). 



(28) 



The additional contribution - £q & u (Q)@o(Q — K) is absent 
in two-component systems. On the other hand, in the normal 
phase T > T c , the SU(3) symmetry is restored and the three 
color becomes degenerate. For the three-component system 
we have 



Z rg (tf) = Eb(K) = -2 2 D (QWo(Q ~ K), 



(29) 



where Do(Q) is the common particle-particle ladder in the 
normal state with A = 0. However, for the two-component 
system, the self-energy in the normal phase is lfl3ll 



E rg (2T) = -YjD (QWo(Q-K). 

Q 



(30) 



Therefore, the self-energy induced by the pairing fluctuation 
effects in the normal phase is enhanced by a factor of 2 in 
comparison with the two-component system. 



which indicates some additional Goldstone modes corre- 
sponding to the broken generators A4, A5, A&, Aj of the SU(3) 
group [24]. They possess a quadratic dispersion u>(q) ~ q 2 at 
low momentum and frequency due to the asymmetry between 
the paired and unpaired colors 12411 . 

Using the above matrix structure of the particle-particle lad- 
der D(Q), we find that the diagonal component of the self- 
energy Zfj (K) takes the form 



Z L U (K) = Z Ig (K)A rg + Z h (K)A h , 



(25) 



where "L lg (K) and 'Lb(K) correspond to the beyond-BCS self- 
energies for the paired and unpaired colors, respectively. 
Their explicit forms are given by 



C. Fermion spectral density 

To study the pairing fluctuation effects on the fermionic 
excitations, we investigate the single-particle spectral density 
function J[{to, k). It can be obtained from the dressed fermion 
Green's function S\\{K). Ensured by the residual SU(2) sym- 
metry, the paired and unpaired colors decouple even though 
the pairing fluctuation effects are taken into account. We have 



Sn(K) = & g (K)A Ig +g h (K)A h , 



(31) 



where Q rg {K) and @b{K) are the dressed propagators for the 
paired and unpaired colors, which can be expressed as 



1 



i(jj n - & -ErgCK^.k) 



ico n +X r;; ( — f'^i ,k) 



Eb(K) = -2YD 5 n (Q)g A (Q-K). 



Q 



(26) 



and 



These results manifest the schematic plot in Fig. [T] Each color 
couples to the other two colors via the corresponding particle- 
particle ladders. The off-diagonal component £p(^0 can be 
evaluated as 



1 



£ k -£b(zw,„k) 



(32) 



(33) 



Q 



(27) 



respectively. We emphasize that the expression for ^(iaj^li) 
holds even though the off-diagonal self-energy I.^ 2 (K) is taken 
into account. In the BCS mean-field approximation the pair- 
ing fluctuation effects are absent and we have @ rg {K) = @ A {K) 
and g h (K) = g (K). 



The spectral functions ^R a (cj, k) for the paired and unpaired 
colors are defined as 

&a((o,k) = — Irn0*(w,k), a = rg,b. (34) 
n 

Here and in the following, the retarded Green's functions are 
denoted by the subscript R, i.e., X R (u>) = X(oj + ie) where 
e — + . In the BCS mean-field approximation we have 



and 



^£ cs (w,k) = u 2 k 6(io - E k ) + ulS(cj + E k ) 



^ cs (w,k) = i5(w-£ k ). 



(35) 



(36) 



Here u\ = (1/2)(1 + &/£ k ) and v\ = (1/2)(1 -&/£ k ) are the 
BCS distribution functions. Since the paired colors obtain an 
large pairing gap ~ A at strong coupling, their spectral density 
function J[ r „{a), k) remains gap-like even though the pairing 
fluctuation effects are taken into account. While there exists 
an additional contribution in the self-energy 2 lg (/T) which is 
absent in the two-component systems, we expect that the qual- 
itative feature of the spectral density ^ rg (w, k) is similar to the 
results for the two-component systems 111311 . In general, pair- 
ing fluctuation and temperature effects lead to broadening of 
the gap-like peaks. 

In the following, we will focus on the spectrum of the 
unpaired blue color, which is unique in the present three- 
component Fermi system. It is interesting to study how 
the pairing fluctuation effects influence the spectrum of the 
unpaired color which takes a free dispersion in the naive 
BCS mean-field description. The spectral density function 
yib(w, k) can be expressed as 



jH b (w,k) = — 



ImEfVw, k) 



:<37) 



n [cj - & - ReE*(w, k)] 2 + [Iml£(w, k)] 

The real and imaginary parts of the retarded self-energy 
Z^(w, k) are related by the dispersion relation 



R f + °° da, 

ReZ*(w,k) = P — 

J-oo X 



+« da/ ImZlW.k) 



(38) 



The imaginary part is explicitly given by 

ImZ*(w,k) = 2^14 _ k ImD??(w + £ q - k ,q) 
q 

x [/(£q_ k ) + Hw + £ q _ k )] 
+ ^_Jm£>^(w-£ q _ k ,q) 

x [/(-£ q - k ) + b(co - £ q _ k )]}, 



(39) 



T=0.05r, 
— k=1.0i„ 



T=r t 

- k=1.0i„ 





where b(E) = 1 /(e E ' T - 1) is the Bose-Einstein distribution. FIG. 3: The spectral density ^ b (w, k) of the unpaired blue color for 



The imaginary part of D x * {a>, q) is given by 



k = k u 



2mji at various temperatures. 



ImD^f (<j, q) = - 



Im^(w,q) 



[Re^R (w ,q)]2 + [ i nw 5R (W;q)] 2' 



(40) 
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FIG. 4: The imaginary part of the retarded self-energy k) of FIG. 5: The real part of the retarded self-energy k) of the un- 

file unpaired blue color for k - k^ = -^2mfi at various temperatures. paired blue color for k = k^ = -\/2mfi at various temperatures. 
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where the pair susceptibility ^[(w, q) reads 



l-/(£k)-/(?q-k) 



/(£k)-/(?,-k) 
i 

£k - £q-k + U 



+■ £q-k - W 
1 

2e k 



(41) 



Now we turn to the numerical results for the spectral density 
J^ip), k). The system is characterized by a single dimension- 
less coupling parameter l/(kpa s ) where kp is defined via the 
total density as n - k^/(2n 2 ). To achieve strong coupling, we 
consider the resonant interaction with a s — > oo (l/a s = 0). In 
this case, all quantities can be expressed via two parameters: 
the chemical potential fi and the effective Fermi momentum 
k M = JTmjx. This also simplifies the numerical procedures. 
In this case, we do not need to solve the number equation © 
to obtain the chemical potential fi in units of the Fermi energy 
Ep = £p/(2m). The order parameter A in units of fi as a func- 
tion of T/fi is solved from the BCS gap equation ( fT2b . Then 
we can calculate the spectral density ^lb(w, k) in units of 1/fi 
through the equations ( |37l i, fl38l l, and ( f39l >. 

The numerical results of the spectral density ^b(<^) for 
fermion momentum k — k^ are shown in Fig. [3] At very 
low temperature T — » 0, we find that there exists a very sharp 
Fermi liquid peak around a> = 0. The spectral weight of the 
continuum part is very small. Therefore, the Fermi liquid 
picture of the unpaired color remains valid even though the 
pairing fluctuation effects are taken into account. This can be 
understood by the behavior of the imaginary part of the self- 
energy T^((d) shown in Fig. [4] Due to the formation of a large 
pairing gap A for the paired colors, the imaginary part van- 
ishes in a wide regime around a> = 0, which can be seen from 
the properties of the Fermi-Dirac and Bose-Einstein distribu- 
tion functions in Eq. (f39|. The real part of the self-energy 
^(w) shown in Fig. [5] leads to a correction to the effective 
Fermi surface, which is shown by the fact that the Fermi liq- 
uid peak is no longer located precisely at to — 0. 

However, as the temperature is increased, the pairing gap 
A is reduced and the contribution from the imaginary part 
of the self-energy Sj^w) becomes more and more important. 
When the temperature is high enough but below the critical 
one T c , the imaginary part ImE^(w) has two main effects: 
(i) It becomes finite in the regime around a> — and hence 
the Fermi liquid peak broadens; (ii) It has two sharp min- 
ima which breaches the Fermi liquid peak and the other two 
gap-like peaks. Therefore, as the temperature is increased, the 
spectral weight of the Fermi liquid peak becomes smaller and 
smaller and there arises two gap-like peaks around the Fermi 
liquid peak. These two gap-like peaks can be referred to as the 
"pseudogap peaks", since they are not induced by the pairing 
fluctuation effects rather than the pairing gap or superfluid or- 
der parameter A. When the temperature reaches T c , the imagi- 
nary part ImS^(w) undergoes a characteristic change: the two 
sharp minima join to a smooth minimum at a> = 0. Therefore, 
the Fermi liquid peak disappears completely at T — T c and 
only the two pseudogap peaks remain. Note that at T = T c the 
SU(3) symmetry is restored and therefore the spectral density 
shown in Fig. [3]is also true for the red and green colors. 



The above discussion shows clearly how the Fermi liquid 
behavior of the unpaired color at low temperature evolves to 
the pseudogap behavior at and above T c due to the pairing 
fluctuation effects. In a wide temperature regime below T c , 
we find that the Fermi liquid peak and the pseudogap peaks 
coexist. It is intuitive to understand this coexistence through 
a naive approximation for the fermion self-energy employed 
in the pseudogap theory of the BCS-BEC crossover 113 , 13111 . 
Since the pair propagator !D 5 U (Q) is peaked around Q = 0, we 
can approximate the fermion self-energy Eb(K) as 



where the pseudogap energy A pg 



K), 
is defined as 



A ps= 2 I>ii (G) - 



(42) 



(43) 



Under this approximation, the propagator @\>(K) can be ana 
lytically evaluated as 

(ico n ) 2 - A 2 



(i(o n - &d[(i<o n Y - £ k 
Therefore, for k = k^, the spectral density M^iaS) reads 



(44) 



A 2 + A 2 



■6{oS) + I 



2 A 2 + A 2 



d(co - ^A 2 + A 2 g ) + 8{io + ^/A 2 + A 2 g ) 



■(45) 



While the broadening effects of the peaks are neglected in this 
naive approximation, this analytical expression clearly shows 
a three-peak structure. The spectral weights of the Fermi liq- 
uid peak and the pseudogap peaks are given by 

A 2 , 1 A 2 g 



A 2 + A 2 g ' 



2 A 2 + A 2 g 



Further, the pseudogap energy can be expressed as 
— fc(w)Im£>^(w,q). 

oo 



(46) 



(47) 



As a result of the Bose-Einstein distribution function b(a>), 
we find that A pg — > for T — > 0. In general, the pseudogap 
energy increases with the increased temperature. Therefore, 
from Eq. [46] we see clearly that the spectral weight of the 
pseudogap peak can be neglected at low temperature and the 
Fermi liquid peak disappears at T — T c where A vanishes. 

Finally, we present a brief discussion about the case that 
the SU(3) symmetry is explicitly broken. Let us consider the 
case that the attractions among the three colors are not equal, 
i.e., gi3 = g23 = g' < g\2 = g- In this case, the pair- 
ing of red and green color is favored. Since the pair excita- 
tions corresponding to the particle-particle ladders D 5 (Q) and 
23 7 (2) (they are still degenerate since gu = #23) are no longer 
gapless, the pairing fluctuation effects become weaker and 
weaker as the coupling strength g' decreases. For g' — > 0, the 
particle-particle ladders D 5 {Q) and D 7 (Q) vanish and hence 
the pairing-fluctuation induced self-energy Xb(G) vanishes. In 
this case, the Fermi liquid behavior persists at any tempera- 
ture, as we expect. 
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III. COLOR SUPERCONDUCTIVITY IN TWO-FLAVOR 
QUARK MATTER 

In this section we study quark color superconductors. The 
perturbative-QCD calculation for color superconductivity [4, 
5] is applicable for quark chemical potentials /i » Aqcd- 
However, here we are interested in the moderate density 
regime of quark matter where the quark chemical potential 
ju ~ 400MeV, and the perturbative approach is not applica- 
ble. Therefore, we adopt a phenomenological four-fermion 
interaction model of QCD J3, 0] . Moreover, since the strange 
quark degree of freedom is not activated at fi ~ 400MeV, 
we consider only the two light quark flavors (« and d). The 
Lagrangian density of our four-fermion interaction model is 
given by 



L = q (i$ + fiyo) q + £mt, 



(48) 



where the interaction part is modeled by a QCD-motivated 
contact current-current interaction 



-Cn 



-G c ^Jqy^qXqy A^q). 



(49) 



quark matter is also characterized by a three-component or- 
der parameter A = (Ai, A2, A3), where 



Ai = —(q T Ciy 5 T 2 A 2 q), 

Q 

A2 = -^(q T Ciy 5 T 2 A 5 q}, 

A3 = ^{q T Ciy 5 T 2 A 7 q). (51) 



Due to the color SU(3) symmetry of the QCD Lagrangian 
(and hence the QCD-motivated model), the effective poten- 
tial 'V(A) depends only on the combination AA 1 = |Ai| 2 + 
IA2I 2 + IA3I 2 . Therefore, like the atomic color superfluidity, we 
can choose a specific gauge A] = A + 0, A2 = A3 =0 with- 
out loss of generality. In this gauge, only the red and green 
quarks participate the pairing and the red-green diquark pairs 
condense, leaving the blue quarks unpaired. In the following 
calculations we will adopt this gauge. 

The Nambu-Gor'kov basis for the present case can be de- 
fined as *P — (q, Cq T ) J . In the Nambu-Gor'kov representation, 
the quark self-energy 1,(K) and the dressed quark propagator 
satisfy the Dyson's equation 



Here q denotes the quark field containing two flavors and 
three colors, and G c is a phenomenological coupling constant 
which should be in principle determined by the vacuum phe- 
nomenology of QCD. In real QCD, the ultraviolet modes de- 
couple because of asymptotic freedom, but in the present four- 
fermion interaction model we have to add this feature by hand, 
through a UV momentum cutoff A in the quark momentum 
integrals. The model therefore have two parameters: the cou- 
pling constant G c and the three momentum cutoff A. 



A. BCS mean-field theory 

The contact current-current interaction ( |49l includes all 
meson-like and diquark-like interaction channels, which can 
be obtained via the Fierz transformation. Since we are inter- 
ested in the temperature regime T ~ O(lOMeV), it is sufficient 
to consider only the J p = + diquark pairing channel. The 
pairing gap and critical temperature of other diquark pairing 
channels are much smaller than that of the J p = + pairing 
channel l33ll . Therefore, we have 



-Cint = j 2j ( < ? ! '?5i"2-"laC<3' T ) (q T C iy 5 T 2 A !i q) 



(50) 



a=2,5,7 



where C = iyoy 2 is the charge conjugate matrix and T\ (i = 
1,2,3) are the Pauli matrices in the flavor space. The new 
coupling constant G > denotes the attraction in the J p = + 
diquark pairing channel. 

Due to the attraction among the unlike colors, at low tem- 
perature the quark matter is a color superconductor. This 
is due to the fact that some gluons obtain Meissner masses 
(color Meissner effect) through the Anderson-Higgs mech- 
anism ll34ll . The color superconductivity in the two-flavor 



S U (K) S l2 (K) 
S 2l (K) S 22 (K) 

% + liy 
If- ny° 



S U (Z) Z l2 (K) 
Z 2 i(#) XiiiK) 



(52) 



Again, from the Green's function relation we obtaib the gap 
equation 



A = - V Tv[iy 5 T 2 A 2 S l2 (K)] 

1 K 

and the baryon number density 



(53) 



(54) 



In the BCS mean-field theory, the quark self-energy is cho- 
sen as 



Z(K) = I, BCS (K) = 



iy$AT 2 A 2 
iy 5 Ar 2 A 2 



(55) 



Again, we set A to be real without loss of generality. In the 
BCS approximation, the dressed quark propagator can be ex- 
plicitly evaluated as 



^ CS (^) 



3a(K)toA is + Q{)(K)T A h , 
Ta(K)t 2 A 2 , 
y 5 S^ s (-K)y 5 , 
y 5 S^ s (-K)y 5 , 



(56) 



where To is the identity matrix in the flavor space. Here the 
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Green's functions @ A {K), T A (K), and @o{K) are given by 



r A (K) 



ia>„ + £ k 



(iu n ) 2 - (£ k ) 2 
iA 



2 A kT0 + 



(ico„) 2 - (£+) 2 



-A k ro, 



(/ w „) 2 - (E k y- 

-■ — ^ A k7o + 
iu„ - £ 



iA 



(iOJn) 2 - (£+) 2 



1 



ia>„ + £ 



(57) 



where = ^(^) 2 + A 2 are quasiparticle dispersions, = 

|k| ± ji, and A k = (1 + jqj ■ k)/2 (k = k/|k|) are the energy 
projection operators for massless Dirac fermions. The physi- 
cal value of the pairing gap A is determined by the BCS gap 
equation 



1 



= 4 I 



l-2/(£~) l-2/(£+) 



IE: 



2Et 



(58) 



The new feature here is the presence of anti-quark excita- 
tions in this relativistic quark system. They generally have 
a excitation gap > /i. At high density and low temperature, 
they are therefore irrelevant degrees of freedom. Regardless 
of these anti-particle excitations, we find similar conclusion 
in comparison with the atomic color superfluidity: the paired 
quarks obtain an excitation gap A, while the unpaired blue 
quarks are gapless and possess a free dispersion in the BCS 
mean-field description. Therefore, the same problem we ad- 
dressed in Sec. II also appears in the present color supercon- 
ductor. 



B. Quark self-energy beyond BCS 

The following considerations are highly parallel to the stud- 
ies in Sec. II. To establish the quark self-energy beyond the 
BCS mean-field approximation, we first construct the particle- 
particle ladder or the "diquark propagator" D(Q). In the color 
superconducting phase, it takes the same form as ( fT3l ), i.e., 



(59) 



The diquark pair susceptibility x(Q) here is given by 

xf 2 (Q) = -\Yj Tx[T2AMQ - K)TlAhSl2{K)] > 

1 K 

xf 2 (Q) = Xu(-Q), xfiiQ) = Xn(Q)- (60) 

Note that the trace here should be taken simultaneously in the 
color, flavor, and the Dirac spin space. With the diquark prop- 
agator 2){Q), the full quark self-energy is given by 



+ lP cs (K), 



(61) 



where the beyond-mean-field contribution 2 L (/T) reads 
2 n(£) = -YjYj D " {Q)T2A * Sn{Q - K)T2Ah > 



a,b Q 



^n(K) = - 2j 2^12(2)^2^12(2- K)r 2 A h , 

a,b Q 

tf 2 (K) = y 5 X L n (-K)y 5 , ^(K) = y 5 ^ 2 (-K)y 5 . (62) 

Due to the same reason we addressed in Sec. II, we use the 
quark propagator of its BCS form S BCS (K) in evaluating the 
pair susceptibility x(Q) an d the self-energy zZ L (K). The pair- 
ing gap A is determined by the BCS gap equation d58l and the 
Goldstone's theorem holds. 

Due to the same symmetry structure as the atomic system 
in Sec. II, we find again that^(g) and D(Q) are diagonal in 



the adjoint space, i.e.,xt 
a = 2 sector, we have 



^and flgf 



Df.6 ab . For the 



X 2 U {Q) = ^-2YTt\_@a{K)6k{Q-K)], 

G K 

X \ 2 {Q) = -i^nr^mQ-K)] 

K 



(63) 



and 



£>n(0 = 



xU-Q) 



x 2 u(Q)x 2 n (-Q)-\x\ 2 (QW 
x 2 u (Q)x 2 u (-Q)-lx 2 2 (Q)] 2 



(64) 



The a = 5 and a = 7 sectors are degenerate due to the residue 
SU(2) symmetry group, and xl 2 (Q) - x\ 2 (Q) - since the 
blue quarks do not participate pairing. We have 

X a n (Q) = [& n {Q)T l 

= F " X Tr [^a(^o(G -K) + 6o(K)g A (Q - AD065) 
G K 



for a = 5,7. Using the BCS gap equation d58D . we find that 
£) a (0 diverges at Q = for a = 2, 5, 7, corresponding to the 
spontaneous breaking of the color SU(3) symmetry 0211 . 

Then the diagonal component of the self-energy "£ L (K) can 
be evaluated as 



X L n (K) = ^ g (K)r A IS + Z b (K)T A b , 



(66) 



where T, r g(K) and l.b(K) correspond to the self-energies for 
paired and unpaired colors, respectively. They take the same 
form as Eq. ( 126*1 1. i.e., 

Q 

X b (K) = -2^!DP n (Q)g & (Q-K). (67) 



The off-diagonal component reads 



SuW = 2^2(0^(2 - K)t 2 A 2 . 
Q 



(68) 
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Above the critical temperature, T > T c , the off-diagonal 
component T.^ 2 {K) vanishes and the self-energies for paired 
and unpaired colors become degenerate, i.e., 1*g(K) = "Eb(K). 
This reflects the fact that the color SU(3) symmetry gets re- 
stored above the transition temperature. We have I. xl (K) = 
I.q(K)tqAq, where Aq is the identity matrix in the color space 
and 



Zo(tf) = -2^£>o(00o(e-£). 



(69) 



Here Dq(Q) — £)jf' 7 (0 is the common diquark propagator in 
the normal phase. The above choice of the quark self-energy 
Z()(^0 is identical to that adopted by Kitazawa, Koide, Kuni- 
hiro, and Nemoto ll35ll in studying the precursor and pseudo- 
gap phenomena of color superconductivity. 



C. Spectral density of blue quarks 

The dressed quark propagator Su(K) including the pairing 
fluctuation effects can be expressed as 



(70) 



where Q rg {K) and Qh{K) are the propagators for paired and 
unpaired colors, respectively. We focus on the excitation spec- 
trum of the blue quark. Its spectral density function ^lb(w, k) 
is defined as 

& h (aj, k) = -^ImTr [y°@*(a>, k)] , (71) 

where the retarded Green's function @£((0, k) is given by 

1 



£ b V,k) = — 



[^( w ,k)]-'-^( w ,k)' 



(72) 



To obtain the explicit form of ^{co, k), we note that the 
retarded self-energy 2£(w, k) has the following Dirac matrix 
structure 

7%(cj, k) = s°(w, k) r ° - r h (u, k) 7 ■ k (73) 

due to the fact that the quarks are massless. The imaginary 
part of the self-energy can be evaluated as 

ImZ°(w, k) = J] [H~(k, q) + H + (k, q)] , 
q 

ki*Ll(co, k) = J] [H~(k, q) - H + (k, q)] (q - k) • kj(74) 
<i 

where the functions H ± (k, q) are defined as 



H s (k,q) = 



E s - s£ s 

q-k ' % qk 5R 



2E* 



X [K^ + ^- k )+/(^_ k )] 

PS , £S 



2E 



ImD^-E' k ,q) 



q-k 



[Hfo-E^+fi-E^]. 



(75) 



The imaginary part of T) x j (at, q) is given by 



Im£)^(w,q) = — 



Im^R( W; q ) 



[Re^( W ,q)] 2 + [Im^( W ,q)] 



(76) 



Here the explicit form of the pair susceptibility x\ x (w, q) reads 

*n(^q) 

= £ - J] [*(k, Q) + /i(q - k, q)] [1 - (q - k) ■ k] 

k 

+ J] [7 2 (k, q) + 7 2 (q - k, q)] [1 + (q - k) • k], (77) 



where the functions J\(k, q) and 7 2 (k, q) are defined as 



and 



" ^k-^qk-^ 



^— ' E, +t , - scj 2E, 

s=± k 'q-k k 
" ^ £ k-^q-k + ^ 



(78) 



(79) 



At the quark chemical potential fi ~ 400MeV, the anti- 
quark degree of freedom is much less important than the quark 
degree of freedom. There for we decompose the self-energy 
I$(oj, k) as 



k) = E b (w, k)A k + k)A+. 



(80) 



Here the minus and plus signs correspond to quark and anti- 
quark degrees of freedom, respectively. The self-energies 
X b (w, k) reads 



2*(w,k) = 2°(w,k)±2£(w,k). 



(81) 



Using these results, the total spectral density ^b(<^, k) can be 
expressed as a sum of the spectral densities for quark and anti- 
quark, 



&b(aj, k) = 3%{oj, k) + 3^{fa, k), 



where 



^>J,k) = - 



1 



ImE. (a), k) 



In [CD-R- ReE!(w, k)] 2 + [ImZ? (w, k)] 2 



(82) 



(83) 



In the following, we focus on the spectral density ^l^ioj, k) 
for quark excitation. 

In the numerical calculations, we fix the coupling constant 
G by using the zero temperature pairing gap Ao through the 
BCS gap equation d58l ). The momentum cutoff A is chosen 
as A = 650MeV. It is generally thought that the pairing gap 
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FIG. 6: The spectral density of the unpaired blue quark for k = \i = 
400MeV at various temperatures T < T c . The zero temperature pair- 
ing gap is set to be Ao = lOOMeV. 



FIG. 7: The spectral density of the unpaired blue quark for k = jj = 
400MeV at various temperatures T <T C . The zero temperature pair- 
ing gap is set to be Ao = 150MeV. 
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A at n = 400MeV is of order of lOOMeV JlH. Therefore, 
the ratio A/fj. reaches the order O(l0~ l ), which is likely in the 
strong coupling regime. 

In Fig. |6]and Fig. [7] we show the numerical results for the 
spectral density ^l h (co, k) at quark momentum k = n and at 
quark chemical potential fi = 400MeV. The zero temperature 
gaps are set to be lOOMeV and 150MeV in Fig. [6]and Fig. [7] 
respectively. The qualitative behavior of the spectral density is 
very similar to that for atomic color superfluidity. At low tem- 
perature, the Fermi liquid behavior of the blue quark persist 
and spectral weight of the continuum part induced the pairing 
fluctuations is rather small. As the temperature becomes high 
enough, these continuum part becomes two visible gap-like 
peaks, or the so-called pesudogap peaks, which coexist with 
the broadened Fermi liquid peak. At and above the critical 
temperature T c , the Fermi liquid peak disappears completely 
and the spectral density show purely pseudogap behavior. On 
the other hand, for weak coupling or small Ao, the pairing fluc- 
tuation effects are rather weak and the Fermi liquid behavior 
persists for arbitrary temperature, as we expected. 



consider here only the spectral density for quarks, the sum 
rule f du>Jl b (u>, k) = 1 does not hold. Instead, we have 

J da>Jlb(to,li) = 1. As the temperature increases, the two 
peaks become more asymmetric and their spectral weights 
become smaller, that is, the spectral density for anti-quarks 
should be important at high temperatures. For T » T c , the 
double-peak structure should disappear, as we expect. How- 
ever, the validity of the four-fermion interaction model be- 
comes questionable for T s> T c . 

Finally, we mention that the naive analytical argument by 
using the approximation like Eq. (l42l in Sec. II also applies 
to the present ultra-relativistic quark system. Therefore, the 
coexistence of the Fermi liquid behavior and the pseudogap 
behavior is quite generic for both atomic color superfluid stud- 
ied in Sec. II and the quark color superconductor studied in 
this section. 



IV. SUMMARY 
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FIG. 8: The spectral density of the quarks for k = fi = 400MeV at 
temperatures T > T c . The upper panel is for A = lOOMeV, while 
the lower for A = 150MeV. 



In Fig. [8j we show the numerical results of the spectral 
density ^HZ(cj) above the critical temperature T c , For T > T c , 
the color SU(3) symmetry is restored and hence the results 
valid for all three colors. We find that the pseudogap behav- 
ior at the Fermi surface k — fi persists up to the temperature 
T ~ (2 - 3)T C for A ~ lOOMeV and /j = 400MeV. Since we 



In summary, we have investigated the pairing fluctuation 
effects in atomic color superfluids and two-flavor quark su- 
perconductors. The common feature of these systems is that 
the pairing occurs among three colors of fermions. Due to 
the SU(3) symmetry of the Hamiltonian, one color does not 
participate pairing, which is unique in the three-color sys- 
tems. This branch of fermionic excitation is gapless in the 
naive BCS mean-field description. In this paper, we have 
established a pairing-fluctuation theory for color superfluid- 
ity/superconductivity by using a T-matrix approach. The the- 
ory has been used to study how the pairing fluctuation ef- 
fects influence the excitation spectrum of the unpaired color 
in these systems. The main conclusions can be summarized 
as follows: 

(i) At low temperature, the Fermi liquid behavior of the un- 
paired color persists even though the pairing fluctuations are 
taken into account. The reason is that, at low temperature, a 
large pairing gap for the paired colors forms which suppresses 
the pairing fluctuation effects for the unpaired color. 

(ii) As the temperature is increased, the continuum part of the 
spectral density, of which the spectral weight is rather small 
at low temperature, evolves to two pseudogap peaks. Mean- 
while, the Fermi-liquid broadens and it spectral weight gets 
suppressed. At and above the critical temperature, the Fermi- 
liquid peak disappears completely and the all three colors ex- 
hibits pseudogap-like spectra. 

(iii) The coexistence of Fermi liquid behavior and pseudo- 
gap behavior is generic for both atomic color superfluids and 
quark color superconductors. The reason is likely due to the 
same symmetry structures of these systems. 
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